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Abstract 

We determine the large-distance behaviour of the static dipole-dipole potential for a wide 
class of gauge theories on nonperturbative grounds, exploiting only general properties of the 
theory. In the case of QCD, we recover the known results in the regime of small dipole sizes, 
and discuss recent nonperturbative calculations. Moreover, we discuss the case of pure-gauge 
theories, and compare our prediction with the available lattice results. 


1 Introduction 

The potential between two static colourless dipoles is the simplest example of interaction be¬ 
tween colour-neutral objects that can be studied in the framework of non-Abelian gauge held 
theories. The main physical application of this quantity is in the study of the interaction be¬ 
tween quarkonia, i.e., mesons made of heavy quarks, which can be treated as static colourless 
dipoles in a hrst approximation. From the theoretical point of view, the study of the static 
dipole-dipole potential poses a nontrivial challenge. Indeed, as one is typically interested in its 
large-distance behaviour to describe the interaction between quarkonia, the interesting prop¬ 
erties of the potential are mainly affected by the nonperturbative behaviour of the underlying 
theory, namely Quantum Chromodynamics (QCD). Another complication stems from the fact 
that the mathematical objects relevant to the theoretical study of the dipole-dipole potential 
are nonlocal operators, namely Wilson loops. 

Calculations of the static dipole-dipole potential available in the literature [D El El m [5] 
deal with the regime of small dipole sizes in SU{Nc) gauge theories, mostly using perturbative 
techniques. Even in this somewhat simpler framework, the determination of the static potential 
requires a careful treatment of colour interactions within the dipoles, and of the nonlocality of the 
Wilson loop, in order to avoid the apparent divergence of the potential. This requires a partial 
resummation of the perturbative series HE], or equivalently a representation of the static dipoles 
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in terms of a series of local operators, in the spirit of the Operator Product Expansion lasiE]. 
For short inter-dipole distances b, larger than the dipole size r but smaller than the typical 
hadronic scale, r <C 6 < 1 fm, one can reliably apply perturbation theory to obtain an estimate 
of the potential, which behaves as Vdd 1/6^ [H El El m [5] . At large distances, instead, one has 
to supplement the perturbative description of the small dipoles with nonpertubative techniques, 
like the chiral Lagrangians used in Ref. [5]. The leading behaviour for 6 S> 1 fm is related to 
the two-pion threshold, and was found to be of the form Vdd ~ e [HE]- 

In this paper we want to study the static dipole-dipole potential in a purely nonperturbative 
setting, starting from the dehnition in terms of a certain Wilson-loop correlation function, and 
using only general properties of the theory, namely its symmetries and its spectrum, to derive the 
asymptotic large-distance behaviour. The basic idea is to insert a complete set of states between 
the Wilson loops in the relevant correlation function, and relate the large-distance behaviour of 
the potential to the spectrum of the theory. 

There are several motivations behind this work. First of all, the fully general results for 
the dipole-dipole potential derived in this paper provide nontrivial benchmarks for approximate 
nonperturbative approaches to QCD, like the anti-de Sitter/QCD (AdS/QCD) correspondence 
or the Instanton Liquid Model (ILM), and to gauge theories in general. In particular, we confirm 
the previous calculations of Refs. [SE], and provide a fully nonperturbative definition of the 
various numerical factors entering Vdd- We also compare our results to the recent determinations 
of Refs. [H [7], based on AdS/QCD and on the ILM, respectively. Moreover, since our results 
apply to a generic gauge theory (with mass gap), it is possible to obtain information on the 
interaction of colour-neutral states in various theoretically interesting limits, like the isospin 
limit, or the quenched limit, and to establish how sensitive it is to these “deformations” of QCD. 

The plan of the paper is the following. After setting the notation in Section [2l in Section El 
we express the static dipole-dipole potential in terms of a sum over a complete set of states. In 
Section E] we study the behaviour of the potential at asymptotically large distances, focussing 
in particular on pure SU{Nc) gauge theory, and on gauge theories with light fermions (which 
include QCD). Finally, in Section El we draw our conclusions. Most of the technical details are 
reported in the Appendices El IHl O and IHl 

2 Notation 

In this Section we briefly summarise the important points concerning Wilson loops, and con¬ 
cerning the sum over a complete set of states, mainly to set the notation. 

2.1 Wilson-loop operators 

In the functional-integral formalism, the Minkowskian Wilson loop LVm[C] is defined as follows 

y^M[e] = ^trPexp|-iff jf A^(A)dA'^| , (2.1) 

for a general path C, where P denotes path-ordering0 and are (Minkowskian, Hermitian) 
non-Abelian gauge fields, taking values in the A"c-dimensional defining representation of the 

^Larger path-times appear on the left. 
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algebra of the gauge group. The case we have in mind is that of gauge group SU{Nc), but 
our formalism extends immediately to any subgroup of the unitary groups. In the operator 
formalism, the Minkowskian Wilson-loop operator reads BM 

VVm[C] = ^trTPexp|-i5£i^(X)dX^| , (2.2) 

where T denotes time-ordering of the (Hermitian) non-Abelian gauge-field operators A^{X) = 

where H is the Hamiltonian operator. In this paper we will be concerned 
only with rectangular paths. In Minkowski space, we will denote by CmIzm, Rm,T) the paths 
running along the contour of the rectangles TZm {cr, r), 

RMicr^r) = zm + Rmct+ Tumt , cr,TE[-i,i], (2.3) 

where 

UM = (0, l,0j_), Rm = {r\\,0,f±), ZM = {b\\,0,b±). 

Notice that here T does not correspond to the time-extension of the loop, which is \r 
For the corresponding Wilson loops (at zm = 0) we will use the following notation, 

>viJ^(r||,rl) = >VM[CM(0,i?M,r)], wiJ^(r||,r^) = >Vm[Cm( 0, i?M,T)] . (2.5) 

The Euclidean Wilson loop for a general Euclidean path C, denoted by VVe[C] in the functional- 
integral formalism, and by VV^[C] in the operator formalism, is defined exactly as in Eqs. (|2.ip 
and (12.21) , except that the fields and the scalar product are now Euclidean, and time-ordering is 
with respect to Euclidean “time”, which is here the fourth Euclidean coordinate. Explicitly, 

^e[C] = ^tr Pexp Iy AEf,{XE)dXEf,'^ , (2.6) 

in the functional-integral formalism, and 

^E[e\ = ^trTPexp|-i5£iB^(AB)dAE^| , (2.7) 

in the operator formalism, where Aea{Xe) = i)Ao(0, and AeAXe) = 

XE)e-^XEA^ i = 1,2,3 (these operator relations must be understood in the “weak” 
sense, i.e., they hold for matrix elements of the operators). The Euclidean rectangular paths 
analogous to those defined in Eq. (12.3p will be denoted by Ce{ze-,Re-,T), and run along the 
contour of the rectangles TZe{<7,e) in Euclidean space, 

TZe{<^,t) = ZE + Re<^+ Tuet , a,TG[-i,i], (2.8) 

where 

= (l,0j_,0), = (0,fj_,r||) = (0,r), = (0, 6j_, fty) = (0,6). (2.9) 


(2.4) 

instead. 
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( 2 . 10 ) 


For the corresponding Wilson loops (at ze = 0) we will use the following notation, 

wP{r\\,ri_) = WE[CE{^,RE,T)], {r\\,rE) = We[Ce{^.Re,T)]. 

The Euclidean and Minkowskian Wilson loops and can be formally related by 

analytic continuation. Indeed, the gauge fields in the Euclidean and Minkowskian Wilson loop 
appear only in the combinations AE^{XE)dXE^ and A^{X)dX^, respectively, which are formally 
related as follows: 

AEii{XE)dXEn = Afj,{X)dX^\r^^^-i^ . ( 2 - 11 ) 

It then follows that 

wP{r\\,r±) = wS^\-ir\i,rE), (2.12) 

again to be understood in the weak sense. 

At a certain stage of the calculation we will need Euclidean Wilson-loop operators running 
along the same paths Ce{0, Re,T) appearing in Eq. (12.101) . but corresponding to a different 
choice of the Euclidean “time” direction, i.e., obeying a different time-ordering. These operators 
will be denoted by 


'^pir\\,r±) = ^trTiPexp i -ig (f AE^iXE)dXE^, \ , (2.13) 

Ac [ Jce{0,Re,T) ] 

where Ti denotes time-ordering along direction I, i.e. (for T —)• oo), along the “long” side of 
the loops. 

Let us finally notice that time- and path-ordering can be disentangled using the continuous- 
product representation for the path-ordered exponential. In the Minkowskian case, Eq. (j2.2p . 
parameterising the path as A(A), with A G [0,1] and A(0) = A(I), one has 


Wm[C] 


Nc 


lim 

N^OO 

31 


E 


T 



(2.14) 


where jk are group indices, 

wW{k) = l-^j^A^{Xk)p, Xk = X[X[k + l)) , Xfc = ^(i(A: + i)) , (2.15) 

with 1 the group identity, and ^{k) in Eq. (j2.14l) are ordered according to X^. Similar 

representations hold for the Euclidean Wilson loop operators defined in Eqs. (12.71) and (I2.13p . 


2.2 Complete set of states 

The approach followed in this paper to determine the large-distance behaviour of the static 
dipole-dipole potential is based on the insertion of a complete set of states in a certain Wilson- 
loop correlation function. We use the complete set of asymptotic “in” states, characterised by 
their particle content, and by the momenta and third component of the spins of the particles. 
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We define here the setup in full generality, so that the results obtained in this paper can be 
applied to a wide class of gauge theories. 

Let the spectrum of asymptotic states contain ngp different species of stable particles, char¬ 
acterised by their mass m(s) and spin S(s), with s G {1,..., Usp}. The particle content of a state 
is specified by the string a = {Ni, N 2 ,, iVn^p} of th® occupation numbers Ng = Ns(a). For 
the vacuum W = 0 Vs we use the notation a = ^. Particles are labelled by a double index ig, 
taking values in the index space S = {ig | G N, s = 1,... ,nsp}. For a given particle content 
a, indices run over the set 


Sa = {is ^ S \ I < ig < Ng{a ), Ng{a) / 0} ; (2-16) 

the total number of particles is A/’q = Ns{a). The momenta, and the third component of 
the spins, 53 *^, of all the particles in a state are denoted collectively as where for a general 
ACS 

^A = {{Pis,S3i,) \ is C A} . (2.17) 

A state is completely specified by a and and will be denoted as follows, 

\^sj = 1 ; in), (2.18) 


where the right-hand side stands for the “in” state with the appropriate particle content. Such 
a state transforms under translations and Lorentz transformations as the properly (anti)symme- 
trised tensor product of the corresponding one-particle states, and obeys the usual relativistic 
normalisation. For off-shell momenta, we denote by (1a = {(Pi^jSsi^) | G A} the collection of 
four-momenta and spins. The total energy of a state is denoted as E(^}s^), where for any A C S' 


E{^a) ^is = \jpl + 




Finally, completeness is expressed as 




where for any ACS 


jdnA = j H 

is GA 


j3 

d Pis ^ 


(27r)32ej 


^3is— 


In the following we will also use the notation 


{{f{^A))W,b = j fi^A). 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 
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Figure 1: The relevant Euclidean Wilson loops. 


3 Dipole-dipole potential from a sum over states 


The potential Vm between two static dipoles of size fi^ 2 , with centers separated by b, is 
from the correlation function of two rectangular T X |Ti^ 2 | Euclidean Wilson loops, in 
of large T: 


T^oo 


obtained 
the limit 


(3.1) 


('T'\ 

where W) 2 '^ = y^E[Ci, 2 ] for properly chosen paths Ci ^2 [see Eq. (12.61) ]. and {...)e denotes the 
expectation value in the Euclidean fnnctional-integral sense. Without referring to a specihc 
Euclidean coordinate system (for reasons that will be apparent shortly), the Wilson-loop con- 
hguration can be described as follows. The size and the relative orientation of the “short” sides 
of the loops and of the separation between their centers coincide with the size and the relative 
orientation of fq, r 2 , and b. The “long” sides of the two loops have length T, lie in the orthogonal 
complement of the vector subspace determined by rq, fq, and b, and are parallel. In a nutshell, 
our approach to the determination of the large-distance behaviour of the dipole-dipole potential 
consists of going over to the operator formalism and inserting a complete set of states between 
the loops. Before setting up the calculation in full detail, let us briefly discuss the potentially 
confusing issue of Euclidean “time”. 

Usually, the long sides of the loops are taken to be parallel to the Euclidean “time” direction, 
so that the loops describe the evolution of the dipoles over an amount T of Euclidean “time”, 
which is eventually taken to infinity. However, the notion of Euclidean “time” is well defined 
only after setting up the Hamiltonian formulation of the theory, while in the Lagrangian for¬ 
mulation employed in the functional-integral formalism the direction of “time” can be chosen 
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arbitrarily, thanks to the 0(4) invariance of the Euclidean theory. In our approach we exploit 
this arbitrariness, and we use two different choices of “time” at different stages of the calcula¬ 
tion. For our purposes, it is convenient at first to take “time” along the separation between the 
centers of the loops. In this way, after going over to the operator formalism, one can extract the 
large-distance behaviour of the loop-loop correlator in the usual way, by inserting a complete 
set of states between the Wilson-loop operators. Schematically, 

(0| (rill, fi^)|n)(n| (r 2 ||, r2x)|0), (3.2) 

n 

iT'\ -> 

where Wg '^(ripf'x) denotes a Wilson loop centered at the origin, spanning a temporal interval 
of size |r|||, and with a long spatial side of size T [for more details see Eqs. (I2.8I1 - (I2.10I) . and 
below in this Section], and a caret denotes the corresponding operator. Here |n) denotes an 
energy eigenstate. The large-distance behaviour of the correlator is then obtained from the 
contributions of the lightest states to Eq. (j3.2[) . To determine the dipole-dipole potential one 
has to further take the limit T —)• oo, which affects the Wilson-loop matrix elements appearing 
in Eq. (|3.2h . As we show below in Section EH these matrix elements can be related to the 
correlation function (again in the sense of the Euclidean functional integral) of a Wilson loop 
with appropriate Euclidean interpolating fields, corresponding to the particles appearing in the 
state |n). To study these correlation functions, it is useful to go over again to the operator 
formalism, but with choosing now the “time” direction along the long side of the loop: 

(0|>vg’gr||,r_L)|n) —> (wg^gry , fj_) (/.^^(xi))^ = (OjTi |wg^gry, rj_) 0£;i(xi)| |0), 

(3.3) 

where (pEiixi) are the Euclidean interpolating fields, and Ti denotes time-ordering of the opera¬ 
tors along the direction of the long side of the loop. Here the Wilson-loop operator is denoted by 
, to make it explicit that a different time-ordering is used [see Eq. (I2.13p ]: although Wg 
and correspond to the same Euclidean path, they are in effect different operators. From the 
representation Eq. (13.3p we can then establish the relevant properties of the matrix element in 
the large-T limit, and by taking T —)• oo we can finally derive the dipole-dipole potential at large 
distances. We want to stress the fact that it is perfectly legitimate to use different choices for 
the “time” direction to recast the same (functional-integral) correlation function in the operator 
formalism in different ways, in order to study different aspects of said correlation function, as 
long as these choices are used consistently. In our case, different choices for the “time” direction 
are made in the study of different correlation functions, namely the loop-loop [Eq. ()3.2p ] and 
the loop-interpolating fields [Eq. (|3.3p ] correlation functions, so that no inconsistency can arise. 
We also want to remark that the physical, Minkowskian time plays no role in our calculation, 
which, starting from Eq. m, can in principle be performed entirely in Euclidean spaceH 
Let us now return to Eq. (EH and proceed in a more detailed fashion. As we have said above, 
the 0(4) invariance of the Euclidean theory allows us to choose freely the global orientation of 

^The Minkowskian Wilson loops used in Section lTTl liave to be regarded simply as a useful mathematical device: 
those loops have in fact no relation with the physical process of static dipoles evolving over a large physical time. 
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the Wilson-loop configuration. For our purposes, it is convenient to choose Ci ^2 as follows (see 
Fig.d]), 

Cl = Ce{ze, Rei,T) , C 2 = Ce{0, Re2,T) , (3.4) 

where the paths Ce{ze, Re,T) have been defined in Eq. (12.Sp . and 

= (0,fi,2±,ri^2||) = (0,ri,2), ZE = (0,6j_,6||) = (0,6). (3.5) 


The Euclidean 0(4) invariance further allows us to set 6 _l = 0 and 6|| = |6 | > 0 with no loss of 
generality. We can thus work in this coordinate frame, and write Vdd = Vdd{b] riy, ri_L, r 2 ||, r 2 _L), 
with 6 = |6 I, so that 

{wr)E{wP)E 

In the operator formalism, this correlation function reads 


G^^Hb;ri\\,fi±,r2\\,f2±) 


(0|T|>vPwf^| |0) 
(0|wP|0)(0|wf)|0) ’ 


(3.7) 


where W 12 — 1^b[Ci, 2] [see Eq. (j2.7l) ]. Eor loops that do not overlap in the “temporal” direction, 
i.e., for 6 > |ri|| | + |r 2 || |, the T-ordering sign can be omitted, and so one can insert a complete set 
of states between the loops. Since in this paper we are interested in the asymptotic large-distance 
behaviour of the potential, we will restrict to this case, without loss of generality. Exploiting 
“time”-translation invariance, we can write 


( 6 ; rill, rix, r2||, r2±) = 4 ?(^; ^i|| > > ^2||, r2±), (3.8) 

where 


_ (3.9) 

= J dCls^e '’®^'^®“^M(^)(05^;ri||,ri_L)M('^)(f2s^;r2||,f2_L), 

and we have denoted as follows the relevant Wilson-loop matrix elements, 


’"111 ’ ^I-L) ’" 211 ,7'2±) = ((m C[') , fi_L)M {Qsc ; r 2 \\, 7*2^) 






(0|>v4(nhEi-)l^gc) 

(0|w4(r||,fu)|0) 


M(^)( 05 ^;r||,rj_) 


(^Scl>^4(nhFi-)|0) 

(0|wf)(r||,fu)|0) 


(3.10) 

where is computed on the path Ce{0, Re,T) and has been defined in Eq. (|2.10p . 

Notice that for the vacuum state G^sJ ~ quantities and can be treated 

at once by noticing that under Hermitian conjugation 











and so it is straightforward to show that 


(3.12) 


In the remainder of this Section we show how the expression Eq. (13.8p for the Wilson-loop 
correlator exponentiates to the form given in Eq. ()3.6p . with the right T-dependence in the 
large-T limit. The strategy we pursue is the following. We first derive, in Subsection 13.11 a 
Euclidean Lehmann-Symanzik-Zimmermann (LSZ) |10l ITT| representation for the matrix ele¬ 
ments, Eq. (j3.10p . and from this we obtain, in Subsection 13.2p a decomposition of the matrix 
elements in connected components, with each component describing, loosely speaking, the in¬ 
teraction of an isolated subset of particles with the loop. This decomposition allows us to prove 
the exponentiation of Eq. (|3.8p . and finally to establish that the correlator exhibits the correct 
dependence on T, in Subsection 13.31 where the final expression for the dipole-dipole potential is 
also reported H 


3.1 Euclidean LSZ representation for the matrix elements 

The relevant Euclidean matrix elements M ; ry, rj_) are related to the analogous matrix 

elements for the Minkowskian Wilson loop, 


by means of analytic continuation [see Eq. (12.121) ]. 

M ; ry, fj.) = {ns ^; -fry, rj.). 


(3.13) 


(3.14) 


Although the physical quantities entering the dipole-dipole potential are the Euclidean matrix 
elements in order to recast them into a LSZ-like expression it is convenient to work 

initially with . The quantity ' admits in fact a straightforward LSZ reduction llOpillj. 
which can be written in the following compact formjl 

= Lim5^n(ii5^)LjJ^(P|^,Ps^;ry,r_L), 

/ (3 15) 
dXs^ , 

®We notice, incidentally, that the exponentiation of Eq. (13.81) could be formally obtained in a straightforward 
way by means of the moments-cumulant theorem. However, this would tell us nothing about the properties of 
the exponent, so that we could not prove that the correlator has the right T-dependence. 

'^The derivation of Eqs. (13.1511 and (13.1611 follows the usual LSZ procedure, the only nontrivial point being 
the definition of a time-ordered product involving local fields and the nonlocal Wilson-loop operator. This can 
however be easily obtained by using the continuous-product representation of the Wilson loop, Eq. (12.141) . which 
allows us to write 

fe=0 i J 

for a general path C, and for a general set of local fields 


T{WM[C]Wcj)i{xi)\ ^ lim 

I iV—>oo 


E 
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(T) 

where ' and 11 are defined as follows 






(0|wif)(r||,f^)|0) 


(^11 > ^±) rii.es, j) 


M 


(3.16) 


y^M\r\\,r±) 


M 


n(^5c)= n 

is eS'cK 


Here we have introduced some notation that we now explain. H denotes the product of the 
“projectors” on the appropriate particle poles and spin components: for example, for a scalar 
particle of mass m, 7 r^^^(p) = — m?] for a spin-^ fermion, 7 r^ 2 ^(p, S 3 ) = {p — m)us^{p) = 

{p‘^—m'^){p+m)~^Us^{p), and so on. Both and H (may) carry Lorentz indices, appropriately 
contracted in the product, and are first evaluated off-shell; the on-shell limit, denoted with 


Lim^^ 


n 

is ^Sa 


lim , 


(3.17) 


is taken after computing the product. The operators are the appropriate local 

interpolating helds for particles of type s, normalised to have free-field one-particle matrix 
elements, i.e., the renormalisation constants required in the LSZ formulae have been absorbed 
in their dehnition, so that are renormalised fields. We have denoted collectively with 
— i^Sa ’ '^‘ 5 “ ^ four-momenta of the particles, and the temporal and spatial components 
thereof. A similar collective notation, Xs^ = {Xg^, Xs^), has been used for the coordinates 
of the local operators, and for the corresponding integration measure, dXs^- In the following, 
when there is no need to distinguish between temporal and spatial components, we do not write 
them as separate arguments. Moreover, Ps^ • Xs^ = finally, in the second 

line of Eq. (j3.16jl we have used the functional-integral representation for time-ordered vacuum 
expectation values, denoting with (.. .)m the expectation value in the sense of the Minkowskian 
functional integral. 

The next step is to Wick-rotate Lm to Euclidean space. By means of a simple change of 
variables, one shows that 


j dXs^e-^^^'^-^^'^cZ\iX^s^,Xs^Mx±)- (3-18) 


By sending ^ ^ e * 2 we then obtain 

Lm^ (P^Pl > Ps. ; i-Pl , Ps.; r||, fu), 


(3.19) 
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where we have introduced the Euclidean quantity 


L^P{PESc4.,PESc',r\ 

{XESc.,^ESc4:,\r\ 


rx) = I dXEsJ^^^-'^^^-cP{XES^,XES^4,-,r\^,rE), 

^ (0|T{wf)(r||,fj.)a,6S.^E 

(0|>vf^(r||,r-'x)|0) 

{^PHr\\,f±) rii.es. ^^^{xEis,XEis4))^ 

{y^P\n,r±))^ 


(3.20) 


where PeSo, = {PeS ai PeS ai) denotes collectively the Euclidean four-momenta PEisj ^ESa = 
i^E SaJ ^E Sai) the coordinates of the local operators, dXESa the corresponding integration 
measure, and PeSc ' ^ESa — Ylis&SaP^i!>^^^Eis^J.■ The Euclidean Wilson loop has been 

dehned in Eq. (12.7p . and are now the appropriate local functionals of the Euclidean helds. 
In the third line we have made contact with the Euclidean functional-integral formalism. 

Inverting the analytic-continuation relation we fincjl 


LP{PLPs^;n,rE) = 

LP{PLPs^-,-ir\\,f±) 


(e-*5 (-POJ, P5.; r|| , r-1), 

= (-i)^“LP(e-*i(-P|j,P5„;r||,r-l). 


(3.21) 


Summarising, is obtained by first computing Lp'^ for real arguments, then performing 

(T) 

the Wick rotation to obtain (with real arguments), taking the momenta on-shell and 

hnally analytically continuing ry —>■ —iry. However, the on-shell projection and the last analytic 
continuation ry —?> —*ry should not interfere. If it is so, then 


M(^)(H5„;ry,r-l) = Lim5.n(H5j(-i)^“Tf^(e-*5(-P|j,P5.;ry,r-l), 


(3.22) 


(T) 

and we can follow a simpler route: after computing for real arguments, we perform the 
partial Wick rotation (e “*2 i~Psa)^ ^|| i finally take the momenta on-shell. 

Let us finally notice that the correlator in Eq. ()3.20p is a renormalised quantity. In¬ 
deed, we are working with renormalised interpolating fields [see discussion after Eq. (13.171) ]. and 

moreover the Euclidean Wilson loop enters through the combination yVp'^/{yVp'^)E, which 

(T) i-TTri 

is a renormalisation-invariant quantity since Wg ' renormalises multiplicatively mM- As a 
consequence, the matrix elements M , Eq. (I3.10p [as well as the Minkowskian matrix elements 

/rp\ 

MP , Eq. (|3.13l) ]. are renormalised (and renormalisation-invariant) quantities. 


3.2 Cluster decomposition of the matrix elements 

The point in relating M T) with the purely Euclidean quantity L is that the latter admits a 
neat cluster decomposition. Furthermore, as the Euclidean functional integral admits a nonper- 
turbative definition through the lattice discretisation, we can perform the formal manipulations 

®The notation {~Psa) indicates that to obtain the correlator at positive (off-shell) energies one starts 
from negative psis^ ~ ~Pis^ nnd then rotates clockwise in the complex psi, 4 plane. 
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('T’\ 

rather safely. To compute the correlation function , Eq. (I3.20p . we can exploit once again 
the 0(4) invariance of the Euclidean theory, and choose the “time” direction as we please. Eor 
our purposes, it is convenient to now take “time” along direction 1, i.e., the direction of the 
“long” side of the loop. Reverting to the operator formalism with this choice for the “time” 
direction, we write 


{Xe = 


(0|Ti <1 wiP(r||,r-l)05.(^E5j MO) 


(0|>viP(r||,fu)|0) 


OsAXeSo^ ^ II ^%\xEis,XEisi) , 

is ^Sa 


(3.23) 


where has been dehned in Eq. (12.131) . and Ti denotes time-ordering along direction 1. As 

Cg ^ is gauge invariant, we can work in the temporal gauge where the long sides of the loop are 
trivial. With this choice of time-ordering and in this gauge, the Wilson-loop operator can be 
expressed in terms of the following Wilson-line operator, 

We{Re) = P exp i^—ig J ^ dXAEij.{XRE)RE)i^ , We{—Re) = We{Re)^ , (3.24) 

where the time-ordering symbol has been dropped, since only gauge fields at = 0 appear. 
In terms oi We, the Wilson-loop operator reads 

VV^P(l,ru) = ^tr {e^^WE{RE)^e-^^e-^^WE{RE)e^^ , (3.25) 

with H the Hamiltonian operator. Since we are ultimately interested in the limit T —oo, we 
consider only the case when T/2 > \xEisi\ Vis £ Sa- Inserting complete sets of states in the 
appropriate sector of the theory (namely, that transforming as a pair of colour charges in the 
fundamental and complex conjugate representation located at a distance Re), we find 

(0|T, {wiArii.rDOsJXEsj} |0) 

|s2)(s2|i?E;ij), 

S1,S2 i,j 

and moreover 

{Re', 'ij\si){si\RE',ij), 
i,j 

where \RE',ij) = [WE{RE)]ij\0) is the “flux-tube” state created by the Wilson line We{Re)- In 
the limit T —)• oo, the dominant contribution comes from the flux-tube ground state, si = S 2 = 
g = g{RE) (since there is a gap with the first excited state), and we obtain 

CE{XES^',r\\,TE) ^ \\^jP{XESA,r\\,rE) = (ffiTi Iff) • (3.28) 


(3.26) 

(3.27) 
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Consider now the case when the interpolating fields cluster into subsets, well separated from 
each other in the “time” direction. More precisely, given a partition AxiSa) of Sa in K parts, 
AxiSa) = {ak}k=i,...,x, consider the limit 

1 ^Ei', il ^ oo , G dk , G ctk^ ^ k A k . (3.29) 

s' 

By appropriately inserting complete sets of flux-tube states between the subsets of interpolating 
helds, one can show that in this limit the sums over intermediate states are dominated by the 
ground state, and so 

X 

C£;(X£:5^;r||,r_L) —)• ?’||,^±) • (3.30) 

Let us now perform a decomposition in connected components in the usual way, i.e., defining 
recursively, for any T, and for A Q S, 

cp^°'^^(^XEA]rii,f±) =cP{XEA]rii,f±)Y H 

K Ak(A)^{A} aeAxiA) 

^3.31) 

where the sum is over all partitions of A, {^d} is the trivial partition, and for 

one-element sets, so that 

cP{XEsA,rii,r±) = Y E 11 cf ^1,fJ . (3.32) 

K AxiSa) <^&^K{Sa) 

In the limit T —)■ oo one has analogously 

CE{XESA,rii,r±) = Y E 11 Cr“(Xi5,;r,,,rl). (3.33) 

K AxiSa) '^^■^K{Sa) 

In this limit, Ce is translation-invariant along the “time” direction, i.e., direction 1, and so, by 
construction [see Eq. (I3.3ip ]. each connected component is also similarly invariant under 

“time”-translations. Moreover, Eq. (I3.30p shows that in the limit T —)• oo, each connected 
component vanishes when at least one of the interpolating helds is very far from the others in 
the “time” direction. Let us make this discussion explicit by writing 

= CT{ta,Xa) , (3.34) 

where ta = ^ is the average “time”-coordinate of the particles in part a, with Na 

the corresponding number of particles, and Xa denotes collectively all the remaining components 
of the coordinates. Here we have dropped the dependence on the dipole size for simplicity. As 
T -A oo, 

lim CT{ta,Xa) = C{Xa). (3.35) 

T^oo 

We can also say something about how this limit is approached. At hnite T, CT{ta,Xa) is 
essentially constant for |ta| ^ f-, and should not change appreciably as long as jtaj < —/tOcorr, 
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where Ccorr is the so-called “vacuum correlation length” [a Hade] and K is some number of 
order 1 , that depends also on the spread of the “temporal” components of the positions of the 
interpolating fields (which again can be at most a few Ocorr since we are considering a connected 
correlation function), but that is independent of T (when T is large enough and only one short 
edge at a time is relevant to this issue). After a transient region of size approximately 2Kacorr) the 
correlator drops essentially to zero when |ta| > '^ + KOcorr- The conclusion is that CT{TTa, Xa) 
tends to a constant function over the interval G the transient regions in terms of 

shrink as T —)• oo, and the slope of the function there diverges. So CritajXa) x(^)C'(^a)) 
or more precisely 

lim CT{TTa,Xa) = x{Ta)C{Xa) , (3.36) 

T^oo 

with x(ra) the characteristic function of the interval [— 5 , 5 ]. 

Cluster decompositions for and can also be written down, in full analogy with 

('T'\ 

Eqs. (I3.3ip and (|3.32l) . Comparing them with the cluster decomposition of Cg , one finds 


Tf)c°nn(PE„;r||,r-l) = J dAg , ““ (r,,, fl), 

M ry, fj_) = Lima JJ (Pis > S3is)]LP ““(e"*5 (-P°), Pa; ry, fj_). 


(3.37) 


isSa 


Here we have made use of the fact that the on-shell projector is factorised. The connected 
components of are easily obtained using Eq. (j3.12p . 


M(^)““(L!a;ry,r-l) = [M(^)““(L!a; ry,-r^) 


(3.38) 


Einally, a similar decomposition can be carried out for the various quantities in the limit T —)• 00 . 

The “time”-translation invariance of for a certain part, a, in some partition, Ak, 

reflects itself in the appearance of delta functions in = limT’^oo imposing the 

vanishing of the total “temporal” momentum of the particles in o. Eurthermore, as vanishes 
when the “time” separation between the interpolating fields becomes large [see Eq. (j3.30p ]. the 
corresponding integration regions give no contribution to and no further delta functions 

of subsets of “temporal” momenta can appear. Einally, as the analytic continuation required 
to obtain the matrix elements M = liniT’^oo Af does not involve PEisi^ these properties are 
inherited by the connected components = limT’^oo Af which contain one and the 

same delta function of the “temporal” momenta as More precisely, for the physically 

relevant quantity M one can write 

= j dtae*''“‘“Fr(ta,Pa;My,M±), (3.39) 


where Qa = Ylis&aPis 1 ) for a certain function Ft, obtained from Ct through integration over Aa, 
Wick-rotation of the momenta, and on-shell projection (see Appendix ICl for more details). The 
important point is that these steps should not change the way the large-T limit is approached, 
i.e., for large T 


lim 

T—^OO 


FT{TTa,Pa\ri\\,ri±) 


x(Ta) AT““(f2a;riy,ri_L), 


(3.40) 
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for a certain from which it follows 

!• __ 7\/T conn__ _-t* \ 


MTill, fi_L) = 5 (Ei.ea^’h l) J’lll, f^l±) • 

For the other connected matrix element, [see Eq. ()3.38p ]. we similarly have 


(341) 


(3.42) 


3.3 Dipole-dipole potential from the Wilson-loop correlator 

The purpose of the analysis of the previous Subsection is twofold. On the one hand, the cluster 
decomposition allows us to write down explicitly the exponential form of the Wilson-loop corre¬ 
lator, Eq. (I3.8jl . On the other hand, the properties of the connected components in the large-T 
limit imply that the correct T-dependence is obtained. 

Let us start from the exponentiation. The decomposition of the matrix elements into con¬ 
nected components is not yet the full story, since what appears in Eq. ()3.8I) is the product of 
the matrix elements and Substituting the cluster decompositions of and 

M(^) in Eq. (13.811 . one thus obtains a double sum over partitions. Each pair of partitions 
•^KiSa), of Sa, with K and K parts, respectively, i.e., AxiSa) = {ak}k=i,...,K and 

= {^k}k=i,...,K^ uniquely rewritten as a partition Tj{Sa) of Sa with J parts, and 

a set of irreducible pairs of partitions \AKj-,Ap^]{Fj) of the parts Fj € Fj{Sa)- By an irreducible 
pair of partitions we mean that there are no proper subpartitions {afc}fc=i,...,A''. C AxjiFj), and 
/f/ C Ax^{Fj)., such that Ufco), = U^a'^. Checking a few examples should convince the 
reader; a formal proof is given in Appendix The double sum over partitions can therefore be 
rewritten as 


E E E E =E E n (eE E 

K Ak(SA R AASA j Tj(Sc)FeTj{Sc) \ k L [AkAAF)/ 

Working out the consequences of this fact is a straightforward but lengthy exercise in combina¬ 
torics, which is described in detail in Appendix [Bj Here we report only the final result for the 
Wilson-loop correlator, which reads 


G^'^Hb;ripri_L,r 2 \\,r 2 ±) = exp J X] n /V ( ib;ripri±,r 2 \\,f 2 ±) I , (3.44) 

where we have introduced the following quantities, 

Qa\b-,riii,fi±,r2ii,r2±) = ^ (( 11 M('^)““(Ha;ri||,ri^) 

K K [AkAkKSA WaeAKiSc,) 

u (3.45) 

X H M(^)““(H,;r 2 ||,r 2 x)\\ 

a&Aj({Sa) II 
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Recalling Eq. (13.61) . the dipole-dipole potential reads 
Vddib-, rill, rix, r2|| ,r2±) = - ^li^ ^ 

^ X/ n ('aV ^ ^ ^ Ric,^K](Sc,)(^i’'i|h^i-L)^2|p^2±), 

K [AkAkKSc,) 

where [AK,Aj^]{Sa) is a pair of irreducible partitions of Sa, and 
- V[A^,AK]iSc)ib’ '■ill ’ ^1^’ '■211 > ^2±) 


= lim — 
T^oo T 


n " 


(^)““"(b!,;ri||,fi^) n M 

{Sa ) {Sa) 


(T)conn 


{^a;r2\\,r2±) 


(3.46) 


(3.47) 




(T) 

The crucial point is now to show that Qa ' diverges linearly with T. As we have argued in the 
previous Subsection, in the large-T limit each connected component develops a Dirac delta of 
the total “temporal” momenta qa = Ylis^aPisi Qa = Si^eaT'Ai in each part. In Appendix 
[^we show that, due to the irreducibility of the pair of partitions, only K + K — 1 < Ma out 
of the K + K linear combinations of momenta qa and q^ are independent, the only relation of 
linear dependence being 


Qa= Y Qa= Y 
a^AK (Sa) aGA^{Sa) isG^Sa 

In practical terms, this means that in the large-T limit the integral in Eq. (I3.47P is divergent, as 
one of the K + K Dirac deltas of Eqs. (I3.4ip and (|3.42l) has to be evaluated at zero. However, 
this also means that the divergence is linear in T, so that it gets cancelled by the 1/T factor, 
and V^^AkAkMSo) i® finif®- ^ detailed calculation showing this, which makes use of the large-T 
behaviour of the connected matrix elements, Eq. (I3.40p . is reported in Appendix O Here we 
quote only the final result. 


- ^[AkAkMSA '■ill ’ ^1^T2|| , ^2±) 

= I dQs. 

X A4““(IIa;D||Ti±) n -^““(^aT2||T2±), (3.49) 

Q'G^xC'S'ck) (S'ck) 

o o 

^AK,AK]{S.M)^HEi,€S^Pisi) n HEi^eaPisi) n 

aSAKiSa) aSA^{Sa) 

where the symbol o denotes that the product is over all the parts in the partition but one. 

The expressions Eqs. (13.461) and (I3.49P fully encode the static dipole-dipole potential when 
the dipoles do not overlap in the direction of their separation, i.e., for all b and fi ^2 such that 
|61 > |ri • 6| -|- |r 2 • 6|. In the next Section we use them to extract the behaviour of the potential 
at asymptotically large distances. 
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4 Asymptotic behaviour of the potential at large distance 


At this point it is straightforward to derive the large -6 behaviour of the potential. From Eq. (j3.49jl 
we see that the 6 -dependence is contained entirely in the factor However, we still 

have to perform the phase-space integration. Making the change of variables Vbpi^ = qi^, we 
can rewrite Eq. ()3.49l) as follows, 


nil, n±, nil, n±) 




CLG^K i^ot) 


{Sa ) 


where 


n j 

r _ 1 o 


d^Qi. 




\ 


1 + 


% 


bml^ ’ 


denotes the sum over the spins of all particles, and we have denoted 

-^CIa = ®3is) I n ^ 4} , A C S'. 

In the limit of large 6 , we can expand £i^ and the integration measure d^Sa S’S follows. 


n. — 1 + 




26m? ^ 

(«) 


+ 0 ( 6 -^), 

d^Qis 


n (27r)32m(,)(l +0(6-1)) 

Is t Oq; 


^dqs^{l+0{b-^)) 


and moreover we can expand the matrix elements around zero momentum, 

= {(0, sai^) I is G A} + 0(6 2) = 0)4 -|- 0(6 2). 

To leading order we find 

3 A/*q; — (“I” — 1 ) L AT / \ 

^ ’ ^1^’ ’’211, n±), 

where the full 6 -dependence is in the hrst two factors, i® ^ constant, 

A 

is^Sa 2m/ 




[■Ak 


= J dqs^e 


\K+K-1 • 


2 -^“ (27r)21^“-(-^+-^-i) 


'('>) (27r) 

JJ dqi^ie 


)(9i) 




-is ^Sa 


[Ak iAp^]{Sc 


.)(9i)’ 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 
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and the dependence on the size and orientation of the dipoles is contained in A4 


[A.K 5-^/c] ’ 


= E n ““ K; ^i|| , n±) _n ^ " 2 || , r 2 ±) ■ 

{sa} aSy^iflSa) a&Aj^{Sa) 


(4.8) 


Here we are implicitly assuming that the connected matrix elements are finite, nonzero 

quantities at zero momentum. This is expected to be the case for states containing only massive 
particlesl^ Notice that the exponent 7 = [3A/q — {K + K — l)]/2 of the power-law term in 
Eqs. and (14.6p obeys the inequality 7 > Ma > 1 (see the end of Appendix]^, as well as 
7 < (3A/q, — l)/2 since K,K >1, for any (non-vacuum) state. 

The leading behaviour of the potential is determined by the contributions of the 

lightest states with nonzero Wilson-loop matrix elements, with higher-order contributions being 
exponentially suppressed. Since the Wilson-loop operator depends only on the gauge fields, it 
is obviously invariant under any symmetry of the theory acting only on the matter degrees of 
freedom. This implies a selection rule involving the corresponding quantum numbers, which 
have to be the same as those of the vacuum in order for the Wilson-loop matrix element to be 
nonzero. In particular, in the case of QCD the Wilson loop is insensitive to flavour, and so its 
matrix elements can be nonzero only for states carrying no flavour quantum numbers, which 
results in a selection rule for baryon number, electric charge, strangeness, etc., that must all 
vanish. 

For the interesting gauge theories, the lightest particle is typically a spin-zero particle. In¬ 
deed, lattice results for SU{Nc) pure-gauge theory indicate that the lightest “glueball” has 
quantum numbers = 0 ^'*' (see, e.g., Ref. m)- For theories with Nf light fermions, the 
lightest particles are the N'j — 1 (pseudo)Goldstone bosons generated by the spontaneous break¬ 
ing of the (approximate) chiral SU{Nf)L x SU{Nf)ji symmetry (at least if Nf is not too large). 
This is the case for real-world QCD (Nc = 3, Nf = 2), where the lightest states are the pions 
(pseudoscalars). For spin-zero particles it is possible to derive easily further selection rules on 
parity and charge conjugation. As we show in AppendixlHl for a self-conjugate particle with C 
and P phases r]c and r]p, nonzero matrix elements are possible only if 77 c = r/p = 1. 

Let us now discuss in detail a few interesting cases. In QCD, the lightest particles are the 
three pions, and tt^, but due to the selection rules on electric charge and on parity, they 
have vanishing one-particle matrix elements. The lightest state with nonzero matrix element is 
the one containing two followed by the state containing a tt+tt” pair. In both cases there 
is a single irreducible pair of partitions contributing to the potential, namely the pair of trivial 
partitions [^ 1 , ^ 1 ], and so denoting with and 8^+^^- the relevant Sa, we find from Fq. ()4.7p 




(dvr) 2 


■^[ApAi]{S^+^.) 


(dvr) 2 


(4.9) 


° In the presence of massless particles, they are expected to vanish, in order to cancel the divergence in the 
phase-space measure. We have verified this explicitly in the simple case of pure 17(1) gauge theory, i.e., for free 
photons. Notice that, in the case discussed in the present paper, the above-mentioned divergence is only apparent 
and does not require the vanishing of the matrix elements. However, one can easily show that in the case of 
Wilson loops at nonzero angle 9, considered, e.g., in Ref. m, there is indeed a logarithmic divergence unless the 
matrix elements vanish. 
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Due to the very small relative mass difference between the neutral and the charged pions (also 
when electromagnetic effects are neglected), the 7r"’'7r“ contribution is appreciably suppressed 
compared to the vr^vr*^ contribution only for distances well beyond the range of the dipole- 
dipole interaction. Therefore, although strictly speaking it is the vr^vr*^ state that determines 
the asymptotic behaviour of the potential, it is physically more meaningful to treat charged and 
neutral pions on the same footing. We will then consider the limit of exact isospin symmetry, 
and ignore the small mass difference between vr^ and vr^. In this limit the contributions of the 
vr^vr*^ state and of the vr+vr” state are identicallll and so, taking into account the symmetry factor 
1/2 for the 7r*^7r^ state, we have to leading ordeiH 


Vdd{b;n\\, ri±, rail, f 2 ±) ^ V^Ai,Ai]{S.n)i^’ "’ll! ’ D^T 2 ||, ^ 2 ±) 


3 ,_e 

—)■ — Jm^y— 
6^00 2 ^ 


6—)-oo 

—2mTrb 


(dvrft)' 


-M (at .; r,||, r,±) M (al „; r,,,, ftx) . 


(4.10) 


Due to Eq. (I3.42p . the potential is attractive at large distancesl^ For Nj degenerate flavours of 
quarks qi, the relevant states are those with pairs of “pions” nijiTji, where vrjj = qiqj for i ^ j, 
and Nf — 1 pairs Trjvrj with TTj a combination of qiqi (the completely symmetric one is excluded). 
There are Nf{Nf — l)/2 pairs with i ^ j, and the Nf — 1 states with two vrj require a symmetry 
factor 1/2; the net effect is to replace 


3 NfjNf -1) Nf-1 N]-l 

2 2 2 2 

in Eq. (I4.10p . 

The dependence on b and the properties of our result, Eqs. (I4.10p and (j4.1ip . agree with the 
hndings of Refs. HE], which apply in the regime of small dipole sizes. On the other hand, the 
calculations of Refs. HE], via AdS/QCD and in the ILM, respectively, report large-distance 
behaviours of the form [see Ref. [ 6 ]i Eq. (37)] and ~ e “”*®^/62 

[see Ref. [7], Eq. (64)], where Mx = mp^jvfjs with mp the rho mass, and mg 35OMeV0 
Concerning the ILM result, the mass scale is approximately of the right magnitude^ On 

^As the Wilson loop is flavour-blind, the only contributions to the Wilson-loop matrix elements M come 
from the isosinglet components of these states, which differ only by a sign. 

^According to the discussion above, in real QCD the strict asymptotic behaviour is obtained from Eq. (14.101) 
by replacing the factor 3/2 with 1/2, and using the 7 r° mass and 7 r° 7 r° matrix elements. 

® More precisely, this is certainly true for small enough dipole sizes if the matrix elements are nonvanishing 
and analytic in and f±_ at zero. If the matrix elements are continuous and never vanish, then this is true for 
all dipole sizes. Furthermore, notice that Eq. (13.4211 implies that this is true for r 2 || = riy, f/jx = —Dx- As we 
show in Appendix El ;r||,rx) depends on rx only through rj, so this is again true for r 2 || = riy, 

|r2x| = |rix|. 

The fractional powers of h look troublesome, since they cannot appear in our general formulas, Eqs. (HU and 
dMj. Indeed, turning around the inequalitites reported after Eq. (14.81) . we have ( 27 -|-l )/3 < AC < 7 , that cannot 
be satisfied by any integer AC for 7 = ^, § • However, in our opinion these values are due to small mistakes in the 
extraction of the asymptotic behaviour of the potential. Correcting these mistakes we find ~ e“^^^/ 6 , 

and ~ e“™''’^/&, which match the form of one-particle contributions to the potential. 

^^The mass ms corresponds to a scalar glueball state in the ILM [7], which is stable to leading order in 1/Nc, 
but which develops a nonvanishing decay width in higher orders [19] , that turns it into a resonant two-pion state: 
this could explain the “anomalous” power-law correction b~^ to the exponential decay in the corrected expression 
for reported in footnote HOI 
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Figure 2: Lattice determination of the static dipole-dipole potential in quenched QCD (data are 
taken from Ref. [22]). Only statistical errors are shown. The dashed line corresponds to the 
asymptotic behaviour, Eq. (I4.13p . with the numerical prefactor adjusted to fit the data points 
at 6 = 0.3, 0.4 fm. 

the other hand, the mass scale Mx is clearly much larger than the pion threshold. However, the 
AdS/QCD correspondence is expected to hold in the large-A^c and strong-coupling limits, and so 
it is perhaps more appropriate to compare this result to the one we have obtained in pure-gauge 
theory, discussed below. In this case the mass scale Mx is of the right order of magnitude, 
although still quite smaller than the lightest glueball mass at large Nc (which is slightly larger 
than at iVc = 3 [18], see below). In the estimate one should probably use the quenched value 
mp for the rho mass, which, however, does not differ too much from the physical value: using 
the quenched lattice results of Ref. [20] for quenched pion masses m)) below m)) < 400 MeV, one 
has mp ~ 800 4- 900 MeV, resulting in Mx ~ 1-2 4- 1.3 GeV, i.e., about 30% below the lightest 
glueball mass. 

In pure SU{Nc) gauge theory, the lightest state contributing to the potential is the one 
containing a single O^"*" glueball, which we denote by 5'o++. In this case there is obviously a 
single relevant pair of partitions, and so 

^[Ai,Ai](So++) = ^ 

SO that 


Vdd{b; Till, rix, r2||, f2±) ^^iMi](5o++) (^; ^i|| > G±T2||, ^2±) 


'* 0 ++'^ 


fe —>00 


dvr 6 


-M' 






S„++T2||T2± 


(4.13) 


Also in this case the potential is attractiveJHl For Nc = 3, i.e., in quenched QCD, the mass of the 
lightest glueball is mg++ ~ 1.73 GeV [5T], corresponding to an interaction range ~ 0.11 fm, 

so that the asymptotic regime should be reached at distances accessible to lattice calculations. In 
Fig.|2]we compare the functional dependence of Eq. (I4.13p with the numerical results obtained on 

'^^See footnote m In Appendix IdI we show that also AI ; ry, fx) depends on r±_ only through rj. 
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the lattice in Ref. [22]. The potential is determined from Wilson loops of length T = 8 and width 
|n, 2 | = 1 in lattice units, on configurations obtained on a 16'^ lattice at /3 = 6.0, corresponding 
to lattice spacing a ~ 0.1 fm. Lattice results and analytical prediction are compatible, although 
within rather large numerical errors. 

The most important subleading corrections come from the expansion in inverse powers of 
b of the energy, the phase-space measure and the matrix elements, keeping fixed the particle 
content, i.e., for two-pion states in QCD and for the lightest glueball state in pure-gauge theory. 
From Eqs. (14.41) and (14.5p . and since terms linear in the momenta in the expansion of the matrix 
elements give vanishing contributions upon integration, we have that the first subleading term 
is of relative order b~^. 

From Eq. (14.6p we see that for a given particle content, with total number of particles J\fa, the 
leading (in b) contribution comes from the irreducible pair of partitions with maximal K + K, 
which cannot exceed Afa + l- In pure-gauge theory, where states with a nonvanishing one-particle 
matrix element are present, the maximal value is attained, e.g., by the pair of partitions where 
one is trivial (the whole set) and one is maximal (each element is a part). In QCD [and in 
similar theories with (pseudo-)Goldstone bosons] there are no such states, and nonvanishing 
matrix elements are at least of the two-particle type. As a consequence, one has K,K < [A/q,/ 2], 
so that K + K < Ma if ■N'a is even, and K K < Ma — 1 if Ma is odd. The leading contribution 
at the A/o-particle level is thus proportional to 


pure-gauge : 


g-A/'amo++^ 


QCD : 








Na even 




Ma > 1, odd 


(4.14) 


It is worth discussing briefly what happens in the presence of massless particles. In this case 
we expect the matrix elements to vanish as powers of the momenta for small \p\ (see footnote 
[ 6 |). Here we drop the particle indices for simplicity. Eor a multiparticle state containing only 
such massless particles, we expect by symmetry that each of them contributes the same power, 
A, of Ip I, to the small-momentum behaviour of the matrix elements. Rather than rescaling 
the momenta as in Eq. (j4.1j) . we now more conveniently set bp = q. Eor large b we find that 
Vdd h with 7 = 1 -|-2A if one-particle matrix elements are nonzero, and 7 = 3-I-4A if matrix 
elements are nonzero starting from the two-particle level. An explicit calculation shows that 
A = 1 for free photons, resulting in the well-known large-distance behaviour of the dipole-dipole 
electrostatic potentially^] If the same value is assumed for massless pions in the chiral limit, then 
we find 7 = 7, in agreement with Refs. BM- 


5 Conclusions 

In this paper we have derived a general nonperturbative formula for the asymptotic large-distance 
behaviour of the potential between two static colourless dipoles, valid for a wide class of non- 
Abelian gauge theories, and for any dipole size. Our result is based only on the symmetries 
and on the nature of the spectrum of the relevant theories, and is therefore a robust result. In 

'^^Notice however that in this case our derivation of the cluster decomposition fails, since there is no gap in the 
spectrum of intermediate states. 
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particular, calculations involving any kind of approximation have to compare successfully to our 
predictions. 

In the case of QCD, we have found the same dependence on the distance as in the results of 
Refs. [HIS], which are valid in the regime of small dipole sizes. We have also compared our results 
to the recent nonperturbative calculations of Refs. [ 61 E], which make use of the AdS/QCD 
approach and of the Instanton Liquid Model, respectively. In both cases we find qualitative 
agreement with our results (apart from some “anomalies” which remain to be clarified). 

We have also discussed the case of pure SU{Nc) gauge theory, for which, to the best of 
our knowledge, there were so far no estimates, and compared our prediction with the available 
lattice results (for Nc = 3) [22], finding agreement (within the rather limited accuracy of the 
numerical data). 

We conclude by observing that the techniques developed in this paper could be easily gen¬ 
eralised to the case of the correlator of two Euclidean Wilson loops forming a nonzero angle 6, 
which is relevant to the study of soft high-energy scattering and hadronic total cross sections 
(see Ref. m and references therein). 
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A Decomposition of pairs of partitions in irreducible subparti¬ 
tions 

Let S' be a finite discrete set. We call irreducible a pair of partitions Ak{S) = {ak}k=i,...,K and 
•^ii'('S’) = {hfc}fc=i K of 'S’j with K and K parts, respectively, if there are no proper subsets 
Is C Ia = {I,---,K} and Is C = {1 ,... , A} such that UfcgXgOfc = ^keXs^k- irreducible 
pair of partitions of S will be denoted by [AK,Aj^]{S). We prove now the following statement: 

Any pair of partitions Ak{S) and Aj^{S) of a set S can be written uniquely as a pair 

Tj{S ), {[AKj,AK.]{Fj)}j , 

where iFj{S) = is a partition of S in J parts, and [AKj,Aj^.]{Fj) are 

J irreducible pairs of partitions of the disjoint sets Fj, with L)j^^AKj{Fj) = Ak{S) 
and U-^^^Aji^iFj) = Ak{S). 

Here the union of partitions of disjoint sets denotes the union of the corresponding families of 
sets. To prove this statement, notice that for any subset S'! C S', a partition Ak{S) provides a 
natural covering of S'!, defined as 

O^[5i] = uf=dafc|afcn5i/0}. (A.l) 
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The following properties of 0_4 hold: 

1. 5i C 0 ^[Si]; 2 . if 5i C 52 C 5, then 0 ^[Si] C ; 

3. OAlak]=ak; 4. if 5i, 52 C 5, then 0^[5i U Sa] = 0^[Si] U . 

Fixed points F = Oa[F\ of Oa coincide with their covering [self-covering), so they must be 
of the form F = Uk^XpO-k for some Ip C Ia- Consider next 0_^a[^] = 
straightforward to show that Oj^[F] = T if and only if F is self-covering with respect to both 
A and A {biself-covering), i.e., F = yJk&XpCLk = ^k^Xp^k fo^ some Xp C Ia and Xp C Ij^. We 
call the partitions {ak}k&Xp and the induced partitions of F. If the induced partitions 

of F form an irreducible pair, then we say that F is irreducible. By definition, an irreducible 
biself-covering set does not contain proper biself-covering subsets. 

The proof now goes as follows. Since, by property 1, F 0'^J^[ak] Vn G N, and 

since 5 is finite, there must be G N such that F^p^ = 0^[afc] = = 0_^AW{k)\: 

i.e., is biself-covering. We now show that the induced partitions of form an irreducible 
pair, so is irreducible. If not, there would be a biself-covering proper subset F' C Ff^p^, and 
since also F" = would be biself-covering, we can assume without loss of generality that 

o^k F F'. Then, by property 2 in Eq. (IA.2p . 

= O^Alak] F Ojj^[F'] =F'C , (A.3) 

which is absurd. A similar argument shows that if ap F F^f.-^, then Ff^p-^ = and analogously 
F{k) = if F with generated from as described above. Finally, the sets 
are all the irreducible biself-covering subsets of 5: if F' is an irreducible biself-covering set, 
then 3ak F F', and by property 2 C F', which contradicts irreducibility unless = F'. 
Obviously UfcE(fc) = 5, and so the set of the J distinct F(j^fs provides the unique 

partition Fj{S) of 5, such that the induced partitions of Fj, denoted by ApjiFj) and Apj{Fj), 
form irreducible pairs [AKj,Aj^^]{Fj). This completes the proof. 

Obviously, to any pair Fj{S), {[AKj,Ap.]{Fj)}j, with [AKj,Ap.]{Fj) any irreducible pair 
of partitions of Fj, corresponds a unique pair of partitions of 5, i.e., ^^(5') = Uj^jAxjiFj) and 
Ap{S) = Uj^iAp.{Fj). The sum over pairs of partitions of a set 5 can therefore be written 
equivalently as 

EEEE =EE n (eE E )■ (a-^) 

K Ak{S) K ApiS) J Tj{S)FeTj{S) \ K R [AkArKF/ 

Consider now the matrices 

= ^ki{i) : = hi{i) ■> (^- 5 ) 

where l{i) and I{i) associate to each element i G 5 the labels of the parts of Ap and Ap^ that 
contain it. The columns and are not all linearly independent, and satisfy exactly J 
independent relations, 

4, j = l,...,J. (A. 6 ) 

{k\ak&Fj} {Hai&Fj} 
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To see this, define the J linear combinations 

Y^^\h,h)^ Y, Y (^-7) 

{k\ak&Fj} mak^Fj} 

which are immediately seen to be linearly independent, as they have no components in com¬ 
mon. There are therefore at most J relations of linear dependence among columns, of the form 
which in components read 

h(€) = h(i) e Pj ■ (A-8) 

We now show that Uj{i) = is constant over each Fj, from which Eq. ()A.6I1 follows. 

By definition, is constant over any a^, and similarly is constant over any a^, and so 
will be Uj{i). Suppose now that Uj{i) is constant over a subset Q Y Fj. Then Uj{i) is obviously 
constant in the covering of Q provided by Ak, since “completes” the parts already present in 
Q. By the same token, Uj{i) will also be constant in 0_^[Q], and in Since Fj = 

for some and € N, applying this argument repeatedly we prove our statement. 

As a final comment, consider the matrix obtained by adjoining the columns {k\ak G Fj} of 
and {k\aj. G Fj} of A^. From the result above, its rank is Kj Y Kj — 1. Since the rank has 
to be smaller than or equal to the number of rows, i.e., the total number of objects in Fj, Nj, 
we have Kj + Kj — 1 < Nj. 

B Exponentiation 

In this Appendix we discuss in some detail the derivation of the exponential formula, Eq. (|3.44p . 
In the previous Appendix we have shown that each pair of partitions AxiSa), Aj^{Sa) of Sa, 
with K and K parts, respectively, can be uniquely rewritten as a partition Fj{Sa) with J parts 
and a set of irreducible pairs of partitions [AKj,Ax.]{Fj) of the parts Fj G Fj{Sa)- Using 
Eq. (|A.4p . and dropping temporarily the dependencies on b, ri 2 ||; and for simplicity, the 
product of two matrix elements M (Qsa) ^ can be written as 



{nsJM(^\nsJ = Y E n 

EE E 

^[AkAk](F)(^f) > , 

(B.l) 


Fj(SYF&Tj{SY 

K K [AkArKF) 

J 


where 

M[AkAk](F)^^f) ^ n ^ 


^(T)conn(j^_)^ 

(B.2) 


a£AK(F) 

a&AAF) 



and we 

have made us of the fact that, by construction. 




P 

II 

Q.a = ^F ■ 


(B.3) 


aeAKF) a&Yii(F) 


A partition Fj{Sa) of Sa in J parts is fully specified by the J strings aj = {Ngj} of occupation 
numbers Ngj = Ns{aj) (the number of elements of type s in part j), satisfying A 0 
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(while Ngj may be zero for some s,j) and = Ng, hy a reference partition with the 

given occupation numbers, and by n“p permutations Pg G <Sns, one for each of the n"p types 
that are present in Sa- This representation is redundant, with J! pairs {{Ngj},{Pg}) 

corresponding to the same partition, since the labeling of the parts is irrelevant and permutations 
of elements of the same type within a part do not yield a new partition. The sum over partitions 
Pj{Sa) can then be written explicitly as 


E 


J! ^ n n,iVs(a,)! 


E 

{Ps£ 5 jVs} 


X] - X] • • • X] n 

{ctjjj 011^9 ^ 


(B.4) 


Consider now the phase-space integral of M Since the integration measure 

is factorised, we have 


(( n ^AK,Ai,]iF){^FW - n ^^[AK,Ai,](F)i^F)J)^ , (B.5) 

\\FeFj{Sc) '' F&Fj{Sa.) ^ 

and so 

= 1^1^ n ^^[AK,A^\{F)i^F)J)^ 

J Fj{Sc)F&Fj{Sc) { K K [AkAkKF) ) 


Taking into account that particles of the same type are indistinguishable, the sum over per¬ 
mutations in Eq. ()B.4p can be carried out trivially, and after a relabeling of the particles we 
get 


G 


(T) 


n. 




E 


1 

J! 


En 

K}5i=i 




-QiP 


(B.7) 


where 

4^^ = EE E ((4^x..4^](5„)(*^sj)) 

K K [AkMMSP 


(B. 8 ) 


with the sum being over irreducible pairs of partitions only. Summing now over states with 
different particle content, and using standard combinatorics results, we finally obtain 


a^^^( 6 ;ri||,fi_L,r 2 ||,r 2 _L) = exp < ^ 






Q^I\h 


;nii,n±,r’2||,?’2±j 


where we have reinstated the full notation. 


(B.9) 
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C Contributions to the potential 


In this Appendix we compute the contributions to the static dipole-dipole potential, 

defined in Eq. (I3.47p . Recall that 


- ^[AkA£{s^) (^; nil > nil, n±) 


= lim — 
T^oo T 


1 


/' 


= lim I dils^e 
T^oo I 


^aGAK(Sa) aeA^iSc) II . j, 

f dtae*''“*“FT(ta,-Pa;n||,n±) ^ ^ 

0,G^K i^Oi) 

^ n / [-fT(ta,-Pa;r- 2 ||,-r 2 ±)]* , 

(‘S'cK ) 

where we have introduced the quantity 
FT{ta,Pa;r\\,f±) = 


isSa 


Lima tt / dAa C'r(ta, A'a; ry, rj.) 


(C.2) 




where dXa = dtadXa, and Pa ■ Xa = Qata + Pa{Xa,Pa), with Qa = I]i^ga?'Li! and ta = 
with Na the number of elements in a. The explicit form of the remainder 
Pa{Xa, Pa) IS not needed. Here we have dropped the subscript E from Euclidean coordinates 
and momenta for simplicity. In terms of Ft we have [see Eq. (I3.39p ] 


M(^)c°nn(Oa;ri||,fi_L) = j dta Frita, Pa] fi±) 

As discussed in Subsection 13.31 for large T 


lim ET(rra,Pa;?'i||Ti±) = X (n) >1 ny Ti±), 

T->-oo " " 


with 


A4““(fla;ri||,rix) = 
Lima t| 

is&a 


dX^P'^P^'^PP C{Xa;rm,riT) 


(C.3) 

(C.4) 

(C.5) 


JP,4^e-'T(-pO) 

The other connected matrix element, [see Eq. (|3.38l) ]. can be similarly recast as 

M ““(Oa; r 2 ||, f2A_) = J di-a [FT{ia, Pa-,r 2 \\, -f 2 x)] * , (C.6) 

where fa = 7 ^ Si^ea^Li and qs. = Yli^eaPisi- la full analogy with what was done above, we 
have 


[Tr(rTa, Pa] r2\\ , -r2±)] * = X (Ta) M ““'"(lla; ?’2||, n±) , 
Vl““"(L!,;r2||,f2x) = [A4““(0,;r2||,-f2x)]* . 


(C.7) 


26 












To compute it is convenient to change variables and use K + K — 1 < Ma out 

of the K + K linear combinations qa and qa, which we denote collectively with q and q, and 
other ‘iMa — {K + K — 1) linearly independent combinations of the momenta, which we denote 
collectively with P. That only + ^ — 1 < Afa of the qa and qs, are independent follows from 
the results of Appendix O Indeed, in the notation of Appendix q = Api and q = Api, with 
Pi denoting collectively all the i, and the matrix obtained by adjoining the columns of A and 
A has rank A' + ^ — 1. The Jacobian of the change of variables can be chosen to be unity, and 
so we can write the phase-space integration measure as 



= / dP / dq 


dq2'K6 


Qa 


Y^q-a 


(C.8) 


where 


dq= H 

aGAK 


dqa 
2tt ’ 


dq= n 


dqa 
27r ’ 


(C.9) 


and J dP is understood to include also the summation over spin, which plays no role in the 
following. We now set 


Prit, t,q, q, P; b, riy, fi_L, r2||, r2_L) 

= e-’^EiUsA JJ FTita,Pa-,rill,ri±) n [Mi-a,P-a;r2\\,-r2±)Y ■ (C.IO) 

Dropping the dependence on 6, and on size and orientation of the dipoles, Pt behaves as follows 
at large T, 

]^im PT{TT,Tf,q,q,P) = P{q,q,P) TT xira) TT xip)- (C.ll) 

aeAK aeAji 

With this notation, and using the integral representation of the Dirac delta, we can write 


- YakAkKsa 

= ^lim ^ J dP J dq J dq J doj J dt J dt .P) ■ 


(C.12) 


Rescaling now qa,qa Qa/T,Qa/T, ta,ta taTppT, and uj -A Tu, and using the large-T 
behaviour of Pt, we hnd 


-14 


[Ak 


j dP j dqj dqj dtj dtj du}P{0,0,P) x{ta 

aeAx aeAg; 

j dP j dwJ^(0,0,P)x(w)^+^ = j dPP{0,0,P). 


^-iqa{ta-aj) 


(C.13) 
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Changing integration variables back to the original ones, this expression can be recast in the 
following equivalent, but physically more clear form, 


- ^[AkA^KSc) ^1|| > ^2||, r2±) 

= f dP f dq f dqT{q,q,P) 27 r 6 {^^qa) JJ 27 r 6 {qa) JJ 27 r 5 (g;j) 

aGAK{Sa) aeAj^{Sa) 

where M[_Ai^,Ag:]iSa) been defined in Eq. (IB.2I) . and 
^[AK,Aii]{Sc,)(Pp = 


HT.is&s^Pisi) n HT.is&aPisi) n HEis&Pisi) 

aGAKiSa) a£Ajf(Sa) 


(C.15) 


where the symbol o denotes that the product is over all the parts in the partition but one. 


D Selection rules for spin-zero particles 

In this Appendix we derive the selection rule rjp = rjc = 1 for the Wilson-loop matrix element 
corresponding to a state with a single, self-conjugate spin-zero particle. To this end, we first 
notice the following transformation laws for the Wilson-loop operator VV£;(r|| ,r±;u) under charge 
conjugation, C, and parity, P [here we have made explicit also the dependence on the orientation 
of the “long” side, up = (1,0,0,0) = (u, 0)]: 

U{C)yVE{r\\,r±-,u)U{C)^ = VVE(r||,fl;-u) = VV£(-'r||,-fj,; u), ^ 

U{P)VVE{r\\,f±]u)U{P)^ = VVE(r||, -fu; -u) = yVEi-r\\,f±]u). 

Notice that “time” is chosen again in the direction of the spatial separation b between the 
dipoles. Under rotations, R, one has in general U{R)VVe[C]U{RP = WeIRC], with an obvious 
meaning of the notation. In particular, for VV£;(r||, ru; u) and for rotations of tt radians around 
the axes r_L and u, denoted, respectively, by R± and Ru, we have 

U{RE)WE{rpr±\u)U{RE)^ =WE{rpr±;-u) = U{C)WE{rpr±;u)U{C)P 
U{Ru)WEin,fE-u)U{Ru)^ =WE{rii,-f±-,u) = U{P)U{C)WEin,fE-u)U{C)^U{P)P 

For states = \p) containing a single spin-zero particle, the relevant matrix element, 

M (p; ry, r_j_), must be of the form 

Af(p;r||,f_j_) =5{p- u) f{p- r±,p^;r\\,rl) , (D.3) 

for some function /, as a consequence of rotation invariance, and of translation invariance along 
u in the limit T -A oo. For a self-conjugate particle with parities pc and pp, one has moreover. 
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from Eq. (ID.2P 


= VcM{R±p;rii,r±) = ?7cM(p;r||,fj_), 
M{p;rii,r±) = T]cr]pM{-Rup;rii,r±) = 'nc'npM{p;r\i,f±), 


(D.4) 


where Eq. ()D.3I) was also used. The selection rules then follow immediately. 

One can further exploit Lorentz invariance of the Minkowskian Wilson-loop matrix elements 
to prove that for spin-zero particles of mass m the Euclidean matrix elements depend only on 
rii and r? in the limit of vanishing spatial momentum. Indeed, for one-particle states 


MM{p;r\\,r±) = 6 {p- u) Fm{p ■ Rm,R\i) = 5 {p- u) FM{p^r\\ - p- r_L, r| - f|), 
and after the Wick rotation ry ^ 

M(p;r||,r_L) = S (p- u) FM{-ip^r\\ -p- r_L, -r| - rl) ^^S{p- ii) FM{-imr\\,-f"^) 
Eor two-particle states, M{pi,p 2 ]r\\,r_[_), one similarly has 

Mm(pi,P 2 ;?’||,^±) = 5 {{pi +P 2 ) ■ u) Fm{pi ■P 2 ,Pi ■ Rm,P 2 ■ Rm,Pi ■ um,Rm) , 
M{pi,p2-,rn,f±) —)■ S{{pi+p2)-u)FM{m‘^,—imrn,-imrn,0,-f‘^). 

p ->-0 


(D.5) 


(D.6) 


(D.7) 
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